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Abstract 

f^ ' In this paper we derive aggregate separation bounds, named after Davenport-Mahler- 

CN . Mignotte (DMM), on the isolated roots of polynomial systems, specifically on the minimum 

r^ ' distance between any two such roots. The bounds exploit the structure of the system 

^ , and the height of the sparse (or toric) resultant by means of mixed volume, as well as 

recent advances on aggregate root bounds for univariate polynomials, and are applicable to 

arbitrary positive dimensional systems. We improve upon Canny's gap theorem [7] by a 

factor of 0(d"~^), where d bounds the degree of the polynomials, and n is the number of 

variables. One application is to the bitsize of the eigenvalues and eigenvectors of an integer 

r) ' matrix, which also yields a new proof that the problem is polynomial. We also compare 

^^ , against recent lower bounds on the absolute value of the root coordinates by Brownawell 

jyl ' and Yap [5], obtained under the hypothesis there is a 0-dimensional projection. Our bounds 

O . are in general comparable, but exploit sparseness; they are also tighter when bounding the 

value of a positive polynomial over the simplex. For this problem, we also improve upon the 

^\j , bounds in [TB] . Our analysis provides a precise asymptotic upper bound on the number 

^ ' of steps that subdivision-based algorithms perform in order to isolate all real roots of a 

^D , polynomial system. This leads to the first complexity bound of Milne's algorithm [52] in 

in 

•£} ' 1 Introduction 

O 

o 



X 



One of the great challenges in algebraic algorithms is to fully understand the theoretical and 
practical complexity of methods based on exact arithmetic. One goal may be towards hybrid 
symbolic-numeric approaches that exploit both exact and approximate computations. Comput- 
ing all roots, in some representation, of systems of multivariate polynomials is a fundamental 
cd ' question in both symbolic and numeric computation. The complexity analysis and the ac- 

tual runtimes typically depend on separation hounds, i.e. the minimum distance between any 
two, possibly complex, roots of the system. This is particularly true for algorithms based on 
subdivision techniques and, more generally, for any numerical solver seeking to certify its out- 
put. Hence, separation bounds are of great use in areas such as computational geometry and 
geometric modeling. 

Davenport [11] was first to introduce aggregate separation bounds for the real roots of a 
univariate polynomial, which depend on Mahler's measure, e.g. [20]. Mignotte [21] loosened the 
hypothesis on the bounds and extended them to complex roots. 
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As for algebraic systems, a fundamental result is Canny's Gap theorem [7], on the separa- 
tion bound for square 0-dimensional systems, see Th. [TUJ Yap |31j relaxed the 0-dimensional 
requirement by requiring it holds only on the affine part of the variety. A recent lower bound 
on the absolute value of the root coordinates [5] applies to those coordinates for which the vari- 
ety's projection has dimension 0, and does not require the system to be square. For arithmetic 
bounds applied to Nullstellensatz we refer to |18j . 

Basu, Leroy, and Roy [2] and, recently, Jeronimo and Perrucci [16] considered the closely 
related problem of computing a lower bound for the minimum value of a positive polynomial 
over the standard simplex. For this, they compute lower bounds on the roots of polynomial 
system formed by the polynomial and all its partial derivatives. This problem is also treated in 

Separation bounds are important for estimating the complexity of subdivision-based algo- 
rithms for solving polynomial systems, that depend on exclusion/inclusion predicates or root 
counting techniques, e.g. [5U1 fT^ 151 [^ fTC] . 

Our contribution. We derive worst-case (aggregate) separation bounds for the roots of poly- 
nomial systems, which are not necessarily 0-dimensional. The bounds are computed as a func- 
tion of the number of variables, the norm of the polynomials, and a bound on the number of 
roots of well-constrained systems. For the latter we employ mixed volume in order to exploit 
the sparse structure that appears in many applications. Any future better bound can be used 
to improve our results. The main ingredients of our proof are resultants, including bounds on 
their height [28] , 

We extend the known separation bound for single polynomial equations to 0-dimensional 
systems, and call it DMM„, after Davenport-Mahler-Mignotte. This improves upon Canny's Gap 
theorem by 0{d^~^). Our bounds are within a factor of 0(2") from optimal for certain systems, 
which is good for n small (or constant) compared to the other parameters. They are comparable 
to those in [S] on the absolute value of root coordinates, but they are an improvement when 
expressed using mixed volumes. It seems nontrivial to apply sparse elimination theory to the 
approach of [S]. More importantly, our result is extended to positive-dimensional systems, thus 
addressing a problem that has only been examined very recently in [5]. 

We illustrate our bounds on computing the eigenvalues / eigenvectors of an integer matrix, 
and improve upon Canny's bound by a factor exponential in matrix dimension. Thanks to 
mixed volume, we derive a bound polynomial in the logarithm of the input size, hence offering a 
new alternative to Bareiss' result [1] that the problem is of polynomial bit complexity. We also 
bound the minimum of a positive polynomial over the standard simplex and improve upon the 
3 best known bounds [21 [IS]j when the total degree is larger than the number of variables. 

Finally, we upper bound the number of steps for any subdivision based algorithm using a 
real-root counter in a box to isolate the real roots of a system in a given domain. This leads to 
the first complexity bound of Milne's algorithm [22] in IR . This aggregate separation bound is 
also useful in the analysis of the subdivision algorithm based on continued fractions expansion 
|19j for polynomial system solving. 

The polynomial systems in practice have a small number of real roots and all roots, real 
and complex, are well separated; it is challenging to derive an average-case DMM„. Another 
open question is to express the positive-dimensional bound wrt the dimension of the excess 
component. 

Paper structure. We introduce some notation, then Sec. [2] derives and proves the multi- 
variate version of DMM as main Thm. [31 Its near-optimality and comparisons to existing bounds 
are in Sec. [3l which also extends it to positive-dimensional systems. Two applications of our 
bounds are in Sec. (H Sec. Ois devoted to subdivision algorithms. 



Notation. O, resp. Ob, means bit, resp. arithmetic, complexity and Ob, resp. O, means 
we are ignoring logarithmic factors. For a polynomial / € Z[xi, . . . ,Xn], where n > 1, deg(/) 
denotes its total degree, while deg^.. (/) denotes its degree w.r.t. Xi. By C{f) we denote the 
maximum bitsize of the coefficients of / (including a bit for the sign). For a G Q, £ (a) > 1 is the 
maximum bitsize of the numerator and denominator. For simplicity, we assume, for any poly- 
nomial, log(dg(/)) = 0{C{f)). Let sep(/), resp. sep(S), denote the separation bound, i.e. the 
minimum distance between two, possibly complex, roots of polynomial /, resp. system (S). For 
/ = o,dY\i=iix — Zi) S C[x], with ad ^ 0, its Mahler measure is M. (/) := 4|ad| ni=i max{l, \zi\}. 

2 The DMM bound 

The univariate case. Consider a real univariate polynomial A, not necessarily square-free, of 
degree d and its complex roots 7j in ascending magnitude, where j S {1,2,..., d}. The next 
theorem [29] bounds the product of differences of the form |7j — 7j|. It slightly generalizes a 
theorem in [2D] , which in turn generalizes [TT] , see also |17[ [T3] . 

Theorem 1 (DMMi). Let f G C[X], with deg(/) = d and not necessarily square-free. Let Q he 
any set of I couples of indices {i,j), ^ ^ i < j ^ d, and let the distinct non-zero (complex) roots 
of f be < |7i| < I72I < • • • < |7d|. Then 

where fred is the square-free part of f. If f £ Z[x], i < d and C (/) = t, then 

^d/2 ^2dr > ^i/2 22£r > "Q |^. _ ^^.| > ^~d ^-d?-Zr{l+d) > ^-d ^-d^-^dr_ 

{i,j)en 

The second inequality follows from: M (/) < II/II2 < (d+ l)||/||oo < {d + 1)^ 2^, e.g. [20l[3l]. 
In the first inequality we can replace Ai (/) by ||/||2- 

The bound of Thm. [T]has an additional factor of 2 wrt |1H I13j. which is, asymptotically, 
not significant when the polynomial is not square- free or d = 0{t). The current version of the 
theorem has very loose hypotheses and applies to non-squarefree polynomials. 

Roughly, DMMi provides a bound on all distances between consecutive roots of a polynomial. 
This quantity is, asymptotically, almost equal to the separation bound. The interpretation is 
that not all roots of a polynomial can be very close together or, quoting J.H. Davenport, "not 
all [distances between the roots] could he had". 

The multivariate case. This section generalizes DMMi to 0-dimensional polynomial sys- 
tems. Let n > 1 be the number of variables. We use x*^ to denote the monomial x^^ • • -x^", 
with e = (ei,...,e„) G Z". The input is Laurent polynomials /i,...,/„, G Er[x^ , . . . , a;„ ] 
= i^[x, x~^], where ET C C is the coefficient field. Since we can multiply Laurent polynomi- 
als by monomials without affecting their nonzero roots, in the sequel we assume there are no 
negative exponents. Let the polynomials be 

fi = J^Cijx''^'^, l<i<n. (1) 

Let {cj^i, . . . , ai^rrii} C Z'^ be the support of /«; its Newton polytope Qi is the convex hull of the 
support. Let MV(Qi, . . . , Qn) > be the mixed volume of convex polytopes Qi, . . . , Qn C IR". 
Here is Bernstein's bound, known also as BKK bound. 



Theorem 2. For /i, • • • /n £ C[x, x~-^] with Newton polytopes Qi, • • • , Qn, the number of com- 
mon isolated solutions in (C*)", multiplicities counted, does not exceed M\/{Qi, . . ., Qn), inde- 
pendently of the corresponding variety's dimension. 

We consider polynomial system 

(S): /i(x) = /2(x) = ... = /„(x)=0, (2) 

where fi € IR[x^^], which we assume to be 0-dimensional. We are interested in its roots in 
(C*)", which are called toric. We denote by Qq the convex hull of the unit standard simplex. 
Let Mj = l\/IV(Qo; • • • ; Qi-i, Qi+i, ■ ■ ■ , Qn), and #Qi denote the number of lattice points in the 
closed polytope Qi. Wlog, assume dim X]"=o Qi — ^ and dim Yli^i Qi — J ^^^ ^^Y I C {0, . . . ,n} 
with |/| = j. We consider the sparse (or toric) resultant of a system of n + 1 polynomial 
equations in n variables, assuming we have fixed the n + 1 supports. It provides a condition 
on the coefficients for the solvability of the system, and generalizes the classical resultant of n 
homogeneous polynomials, by taking into account the supports of the polynomials. For details, 
see [9]. 

Let D be the number of roots G (C*)" of (S), multiplicities counted, so D < Mq. We also 
use B = {n-1) (^) , and dg(/,) = d^ < d. If /, G Z[x=ti], £ (/i) = r, < r, 1 < i < n. Now vol(-) 
stands for Euclidean volume, and (i^Qi) for the number of lattice points in Q^; the inequality 
connecting {jf^Qi) and polytope volume is in [3]. We present the abbreviations and inequalities 
used throughout the paper: 



D<MQ<J\^di< d", B < nD^ < n JJ d- < nd 

i=l i=l 

n n 

(#Q») < n\ vol(Qi) + n < d^ + n < 2d^ 

n 



an 



(3) 



.,2_ 



i=l 

h<{n+ l)^e <{n+ if"r<'"~U^'''"-\ 

n n 

g = J] (#Q,)A^« < 2^"=! ^» n df' < 2"''"" d^ 



"n—1 ,^ 2 Jn — 1 



i=l i=l 



Theorem 3 (DMMjj). Consider the 


0-dimensional polynomial system (S) in 


(0). Let D be 


the number of complex solutions of the system in 


(C*)'", which are 


< |7i| 


< |72 


< ... < 


7d • Let Q be any set of i couplet 


of indices {i,^ 


j) such that 1 < 


i < j < D and 


7i / 7i- 


Then the following holds 
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n){D^+D{e 
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{i,j)en 










(4) 



where \Ur\ denotes the discriminant of the square-free part of the u— resultant, and \ ■ \ 
denotes absolute value. If fi G Z[x] and jj^k stands for the k-th coordinate, 1 < k < n, of 
'-fj, then: 

(2^^C)-^<|7i,fe|<2^^C, (5) 

sep(S) > 2-(3D+2)(D-i)/2 (/dTT^C)-^. (6) 



The following corollary employs mixed volumes. 
Corollary 4. Under the hypothesis of Th. O for fi G Z[x^^], i = 1, . . . ,n, we have 

2M0 (l+Mo+Er=i M, (r+lg(#Q,))) 

2-2 Mo Er=i M, (T+lg(#Q,)) -2 Mg (1+ lg(n+l)+n lg(n)+2n Ig(Mo)) ^ 
2-(Mo+Er=iM.{T+lg(#Q,))) < l^.^l < 2^0+EILiM,{r+lg{#Q0) (7) 

sep(S) > 2-^o(i^»+^§(^")+^'=i^'(^+^§(#Q'))). (8) 

Corollary 5. Under the hypothesis of Th. O for fi G Z[x=''-'^], clg(/j) < d and C (fi) < t, we 
have 

TT 1^. _ ^.| > 2-(3+41gn+4nlgd)d2" 2-2n.(l+n.lgd+r)d2"-l^ ^g^ 

(ij)eo 

2-d"-n(r+nlgd+l)rf"-l < 11 <; 2'i"+"(^+"lS<^+l)'i""^ (10) 

sep(S) > 2-2rf'"-"(2"ig'^+-)'^'"-\ (11) 

Proof of main theorem. Let us first establish the lower bound. Let 7j = (7j,i, . . . , 7j,n) G 
(C*)", 1 < i < D, be the solutions of (S), where fi are defined in ([1]), We denote the set of 
solutions as y C (C*)"'. We add an equation to (S) to obtain: 

(So) : /o(x) = /i(x) = • • • = /,(x) = 0, (12) 

where 

/o = u + rixi + r2X2 H h r^x^, (13) 

ri, . . . , r„ G Z to be defined in the sequel, and u is a new parameter. Now u = — ^^ 7'j 7j^j, on 
a solution 7^-. We choose properly the coefficients of /o to ensure that the function 

/o : y ^ C* : 7 ^ /o(7) 

is injective. The separating element shall ensure injectivity [31 El [TU [27] . 

Proposition 6. Let V C C" with cardinality D. The set of linear forms 

T={ui = xi+ix2 + ---+i''-^Xn\0<i< B={n-l)(j} 

contains at least one separating element, which takes distinct values on V. 



Corollary 7. For /o G -F it holds that ||/o||oo < ^"~ , and 



fo||oo<||/o||2<2i?"^i = 2(n-ir-i 



D 



n-l 



Proof: The first inequality is evident from tlie definition of infinite norm. For tlie second 
inequality, B = {n — 1) (^2)- 



11/. 



Oiloo 



< ll/olb < VT+WTB^+^^^+Wy^ 
<./^<J^^<Vab^^ = 2B--K 



D 



We consider the u— resultant U of (Sq) that eliminates x. It is univariate in u, with coeffi- 
cients homogeneous polynomials in the coefficients of (Sq), e.g. [9]: 



,k n-k Ml M2 



,M, 



U{u) = --- + QkU^r^ '^^Tk^'il ■ ■ ■ <l + • • • 



(14) 



where Qk G Z, c-^ denotes a monomial in coefficients of fj with total degree Mj, and r^ ~ 
denotes a monomial in the coefficients of /o of total degree D — k. The degree of U, with respect 
to n is D. It holds that 



f,Mi Ma M„ 



<c=nii/» 



Mi 

«lloo ' 



(15) 



j=i 



Prom Cor. [7] we have that jr^l < ||/o||oo ^ -B" ^, for all A;. Let |^fc| < /i, for all k. Then 
using [25], see also Eq. (J3D, we get that 

n n 



i=0 



We can bound the norm of U: 



i=l 



i"^ 112 



„ ^D~k^Mi^M2 






D 

< E 

fc=0 

< ^ |/l(S"-l)^-'^C7|^ < /l^C^ ^ (_B2n-2)D-fc 
fc=0 



E 

fc=0 



< /.^ C2 f; (i?2-2). < ^,2 ^2 (^^-^^)^+;-l 



fe=0 



and so 



< eC^A ^B^n~2^D < 4 /j2 (^2 52(„-l)D^ 

|t/||oo < ||f/||2 <2/iC5(""i)^ < 2(?i + l)^^CS("^i)^. 



If Uj are the distinct roots of U, then by recalling the injective nature of /o, we deduce that 
Uj = — Er=i '^« 7i,*- Actually the u— resultant is even stronger, since the multiplicities of its 
roots correspond to the multiplicities of the solutions of the system, but we will not exploit this 
further. 



Proposition 8 (Cauchy-Bunyakovsky-Schwartz). Let ai, 02, . . . , a„ G C and bi,b2, ■ ■ ■ ,bn G 
C Then, 

|ai 61 + • • • + a„ bnf < (|aip + • • • + |a„p) (|6ip + • • • + |6„p) , 

where en denotes the complex conjugate of ai, and 1 < i < n. Equality holds if, for all i, ai = 
or if there is a scalar A such that bi = Xa^. 



Consider 7j,7j and let Ui,Uj be the corresponding roots of U. Using Prop. El 

|n(7J,i - 7j,i) -^ ^ rn{-fi,n - 7i,n)|^ < (r? H h r^) (|7i,i - 7j,ip H h |7J,n - 7j,nP) <^ 



Yl n li,k -"^rk 7j 

k=l k=l 

and thus 



- Y^k-Yhi,k-ij,k\^ ^\'^i-uj\^ ^ Z^^fe ■l^»~^jl^' 



fc=l fc=l \fc=l 

/ n \ -1/2 



^fc=l 



To prove the lower bound of Th. [3l we apply the previous inequality for all pairs in ft, 
\fl\ = i. So we get 



n \ 2 



n i7*-7ji> [Y.^'k] n i^^-^ii- (16) 

(j,j)eo \fc=i / {i,j)&n 

It remains to bound the two factors of RHS of the previous inequality. To bound the first 
we use Cor. [3 It holds ELi ^fc < 1 + T.k=i4 < ll/olll < 4-B2"-2, so 

n 

(^r2)-|^>2-^B(i-")^ (17) 

fc=i 

For the second factor of ()16p we apply DMMi to f/; and thus 



rr iif 7/-I > 9^--D(£'-i)/2 iir/iii-i'-^ /rfn 

Combining (J16p with (J17p and (jlSp . we have the lower bound. In the case where the polynomials 
are in Z[x], then it holds that the absolute value of the discriminant of a square- free polynomial 
is > 1, and we can omit it from the inequality. If the polynomials are in Q[x] the bounds are 
almost the same, since they depend on Mahler's measure. 

Let us now establish the upper bound. We specialize /o in (fT3]) by setting r^ = — 1, for some 
i £ {1, . . . ,n}, and rj = 0, where 1 < j < n and j ^ i. Wlog assume ri = —1. We compute the 
ti— resultant of the system, which we call TZi E Z[u]. Its roots are the first coordinates of the 
isolated zeros of the system, viz. 7i^j, 1 < i < D. Thus dg(7^i) < D. 

The coefficients of 7?.i are of the form, gkC^^C2^ ■ ■ ■ c^^" , where g^ £ Z and the interpretation 
of the rest of the formula is the same as in the previous section. Using |28j . see also Eq. (JH), we 
get that 

n n 

since now /o is a simplex in dimension 1. It also holds that jc^^ ^C2 ^ . . . cj^"| < C. Combining 
the two inequalities we deduce that 

ll^illoo <2^gC, 

and also ||7^l||oo < ||■7^l||2 < 2^ VD + lgC. 

From Cauchy's bound for the roots of univariate polynomials, e.g. [20], we know that for 
ah the roots of 7^l it holds that {2^ gC)'^ < l/||7^l||oo < |7ij| < ll^ilU < 2^ gC. The 
inequality holds for all the indices i and j. Hence, all roots of the system in (C*)" are contained 



in a high-dimensional annulus in C", defined as the difference of the volumes of two spheres 
centered at the origin, with radii 2 gC and (2 gC)"^, resp. This proves Eq. ([5]). 

Now we are ready to prove the upper bound of Eq. ^ in Th. [31 For all a, 6 € C it holds 
that 

|a-6| < 2max{|a|,|6|}. (19) 

Let the multiset Q = {j \ {i,j) G il}, where |0| = i, then 

n l7^ -7.1 < 2^ n l^^-l ^ 2^(2^£'G)' < (2^^^^C7/. 
(«j)en j(zQ 

For proving ^, let {i,j) be the pair of indices where the separation bound of (S) is attained. 
Then 



sep(S) = |7i-7il 



. ^ {li,k - lj,kY > hi,i - 7i,i| > sep(7^l), 
\ fc=l 



where k is any index such that 7^^ / 7^-^ and sep(7^i) is the separation bound of TZi. An 
easy bound on the latter can be derived by applying Th. [1] to TZi with £ = 1: sep(7^i) > 

2l-(?)||7^l||2^ > 2(i-^)(^+2)/2(Z) + l)-^/2(2^^C')-^ > 2-(3^+2)(^-i)/2(/DTTeC7)-^, 
which completes the proof of ([6]) . 

Remark 9. It is tempting to try to prove the lower hound of Th. [3| by applying DMMi to TZi, 
instead of U, as we did in the previous section. This would allow us to eliminate the factor 
j^-(n-i)(D +i(D+i)-D) fj-Qjj^ ijjQ result. However, if we apply DMMi to TZi, it is not obvious 
that the requirements of Th. [T] hold, i.e. that the ordering of (the coordinates of) the roots 
is preserved. Moreover, the bounds on the u— resultant are of independent interest, since the 
latter is used in many algorithms for system solving, e.g. {3 [131 E^- 

3 Comparisons and extensions 

One of the first multivariate separation bounds was due to Canny, later generalized to the case 
when only the affine part of the variety is 0-dimensional [31] . 



Theorem 10 (Gap theorem). -^^t /i(x), . . . , /n(x) be polynomials of degree d and coeffi- 
cient magnitude c, with finitely-many common solutions when homogenized. If^j G C" is such 
a solution, then for any k, either ^j^^ = or \^j^k\ > (3dc)~" 

Let C {fi) = r, then this becomes 2(^S'^+'s'^+'^)"'^", which is worse than the bound in Eq. (jlOp . 
by a factor of 0((i"~^). Li [5], they only require that the system has a 0-dimensional projection; 
m is the number of polynomials and b < n the dimension of the prime component where the 
0-dimensional projection is considered. The bound is: 

l7.j| > ((n + l)2e"+2)-"("+i)'^"(6"-^-im2^)-("-'')'^"'''"\ 

This is similar to ours in ([5|), and we make a comparison in the sequel. Moreover, Cor. [5] 
does not depend on the (total) degree of the equations, but rather on mixed volume, which is 
advantageous for sparse systems. 

A natural question is how close are the bounds to optimum. Let us consider the following 
system [7]: 

2'^x'l = xi, Xj = Xj_i, 2 < j < n. 



The roots are Xj = {2~^'f \ for 2^ » 1. Th. [3] implies Xj > 2-'^"-"(^+"is'^+i)'^" \ which, 
if r ^ (i, is off only by a factor of 2" asymptotically. The negative exponent of our bound 
is 0{n{nlgd + T)d"~^), Canny's bound gives a negative exponent of (Ig 3 + lgd + T)n(F' = 
0{nT(F'). The bound from [5] has negative exponent: n{n+ l)(21g(n + 1) + n + 2)d" +n(lgn + 

We now consider the case that (S) is not 0-dimensional. Then, the bounds of Th. [3] do 
not hold because they are based on bounding the infinite norm of the n— resultant, which is 
identically zero. Specifically, the (sparse) resultant vanishes identically when the specialized 
coefficients of the polynomials are not generic enough, i.e. the variety has positive dimension, 
or, simply, if the variety has a component of positive dimension at infinity, known as excess 
component. 

To overcome the latter. Canny introduced the Generalized Characteristic Polynomial (GCP) 
[8] for dense systems. We use its generalization, called Toric GCP (TGCP) [10]. We consider 
(So) in P^ and perturb it: 



(So) 




0, 1 < i < n. 



where pi = Xlaei" s'^'^'^^x*^, uJi{-) are (suitable) linear forms, s a new parameter, and Pj is 
the subset of vertices in Qi corresponding to monomials of fi on the diagonal of some sparse 
resultant matrix; at worst, Pj contains the vertices of Qi. This perturbation does not alter the 
support of the polynomials nor the mixed volume of the system. 

The TGCP is the sparse resultant of (Sq), denoted T e (Z[c,r])[ti, s], where c corresponds 
to the coefficients of fi and r to the coefficients of /q. The lowest-degree nonzero coefficient of T, 
seen as univariate polynomial in s, is a projection operator: it vanishes on the projection of any 
0-dimensional component of the algebraic set defined by (Sq). We call this Tu G Z([c,r])[u], 
and dg{Tu) < Mq. The roots of Tjj are the isolated points of the variety plus some points 
embedded in its positive-dimensional components. It remains to bound the coefficients of Tu. 
Repeating the construction of U in Eq. ([HI), we get 



where pk € Z, and c^ ^ is a monomial in the coefficients Cij,s, of total degree Mj. It is an 
overestimation, wrt the height of T, if we suppose that Cj ^ is obtained by adding s^ to each 
coefficient of Cj^^, where A = maxj^a{'^i(ci)}. If we expand c-^, the absolute value of the 
coefficients of s is bounded by (^^^2) ll/»ll^' ^ '^^'WfiW^' /V^l. If we expand the term tf, of T, 
the degree of s is bounded by A • HiLi '^«i ^^'^ ^^^ coefficients are bounded by 

n n 

n 



\Qk\ 



n M^ • l^.l • Ir.r^-'^ • n 2''1I/^II^7^/M; = 

n 



i=l 



since every factor c^^, contributes at most Mj coefficients. The bound holds for (the absolute 
of) all the coefficients of T if we consider it as bivariate polynomial in s, u. Recall that \Qk\ < h, 
for all k, where h is defined in Eq. (j3|). This expression also defines A, C. 



Now k < Mq. If we consider Tjj as a univariate polynomial in s, then its coefficients are 
univariate polynomials in u, with degree < Mq. For the 2-norm of Tjj, we use a summation as 
in the 0-dimensional case, and get 

\\Tu\\o.<\\Tuh<2hACB^^-'^'^°. 

The previous bound is the one on U multiplied by A. Thus we can provide a theorem 
extending Th. [3]to positive-dimensional systems, by replacing C by AC, in Th. O 



Theorem 11 (DMMjj with excess components). Consider the polynomial system (S) in 
(0), which is not necessarily 0-dimensional, and where it holds that fi G Z[x], dg(/j) < d, 
and C {fi) < T. Let D be the number of the isolated points of the solution set in (C*)", 
which are < |7i| < I72I < • • • < \jd\- Let fi be any set of i couples of indices {i,j) such 
that 1 < i < j < D, and jj^k stands for the k-th coordinate of-jj. Then the following holds 

{2^'>+^ gCAY > Y\ hi-lj\ >2"^"(^""^)('^"+2)/2(/iCA)^"^"~^S(^^")(^«+^"(^^^)+^\ 
(i,j)en 

{2^° gCA)-^ <\jj^k\<2^''gCA, (20) 

sep(S) > 2-(3Mo+2){Mo-i)/2 (^M^^^^c A)-^«, (21) 

We also have the following, less accurate bounds: 

{i,j)en 

sep(E) > 2-("''-")'^"is^-2'^''"-"(2"ig'i+^+i)'i'""\ (24) 



4 Applications 

We illustrate the bounds of Th.[3]in two applications. The first concerns matrix eigenvalues and 
eigenvectors, and is a standard illustration of the superiority of mixed volumes against Bezout's 
bound. The second is lower bounds of positive multivariate polynomials, inspired by [2]. 

Eigenvalues and eigenvectors. Consider an n x n integer matrix A, with elements < 2"^. 
We are interested in its eigenvalues A, and its eigenvectors v = (fi, . . . , Vn)~^ ■ This is equivalent 
to solving fj = Yll=i ^ij'^j ~ -^^i' ^^i^n, l<j<n, and fn+i = Y17=i '^f ~ ^- ^^ have 
||/j||oo ^ 2"^, ||/n+i||oo ^ 2. The Bezout bound is 2""*"^, whereas the actual number of (complex) 
solutions is 2n, which equals the mixed volume, e.g. |14] . 

Canny's Gap theorem (7j implies \z\ > (6 • 2'^)"''"'+^)^ , for any eigenvalue or eigenvector 
element z ^ 0. Thus, we need 0(nr 2") bits. We get the same exponential behavior in n if we 
apply [31] or [5]. 

It is reasonable to assume that the system is 0-dimensional and apply ([5]) of Th. [3l It 
holds that Mj = 2n, M^+i = n, (#Qn+i) < 2"+^, and {#Qi) < 2"+^ where 1 < j < 



10 



n, and C = ||/n+i||^"+^ W^iWfjC' < 2^^W^ 2" = 2^^'-+-, g < U1=i i#Q^)''' < 
(#Q„+i)^"+i nr=i {#Qi)'''; hence g < (2^+2)" nr=i(2"+')'" < 22-^+5-^+2-. 
The solutions lie in C"+^. The lower bound of Th. [3] yields 

11 ^ 9-2n3-5n2-5-2n2r 

\ \ '^ 1 

where z is an eigenvalues or element of eigenvector. This is exponentially better than the 
previous bounds. Eq. ([6]) from Th. [3] bounds the system's separation bound: — lg(sep(S)) < 
4n^T + nlgn + 4n^ + lOn^ + 12n2 + n— 1 = 0{n'^ + n'^T). This is polynomial in the size of the 
input, and hence we obtain a new proof of Bareiss' result [1], that computing the eigenvalues 
and eigenvectors of an integer matrix is a polynomial problem. 

Positive multivariate polynomials. We consider the following problem, studied in [2]. 
Let P G Z[xi, . . . ,x„] be a multivariate polynomial of degree d which on the n-dimensional 
simplex takes only positive values. We are interested in computing a bound on its minimum 
value, m. We may assume that the minimum is attained inside the simplex; if not, apply 
a transformation which slightly changes the bitsize of P [2]. Let r bound the bitsize of the 
coefficients of P. We wish to find compute a lower bound on tti, greater than zero, depending 
on n,d,T. Equivalently, we have a system with unknowns m,Xi: 

■g^{xi,. . . ,Xn) = • • • = g^(xi,. . . ,Xn) = 0, ,^^, 

P{xi,... ,Xn) = m. 

We use Th. [TTJ since there is no guarantee that the system is 0-dimensional. However, Th. [TT] 
provides bounds for the isolated points of the variety. Since the minimum could be attained on 
a non-zero dimensional component, we should argue that the bounds take care of this case. We 
consider all the irreducible components of the variety defined by ()25p . Each of them contains 
a point for which the bounds of Th. [11] apply. Such a point is the limit of a solution of the 
perturbed system depending on the parameter s when s — ?• 0. Moreover, it is a zero of the first 
non-zero coefficient Tu, seen as a polynomial in s [8l[T0]; Th. [11] bounds these zeros. Now, on 
each of these components, the value of m is constant, since the gradient of P is 0, and so the 
bounds apply for it as well. 

Let Pi = §§- and P„+i = P - m. It holds that deg(P„+i) = d, deg(Pi) <d-l, \\Pn+i\\oo < 
2^ llPilloo < dllVn+illoo < d2^, M„+i <{d- 1)", Mi < d{d - 1)— 1, and D < Mq < did - 1)". 
Using (|20|) we deduce 1/m < 2^ gC A. It remains to bound the various quantities involved, 
defined in (JH): 

n+l n+l 

r* ^ TT II pw^i — II p , i|M„+i TT II p.||Mi 

- ii ll-^«lloo ~ iKn+llloo II IP illoo 

i=l i=l 

n 

\ri— 1 



^ 2{n+l)Td{d-l)"-^ +nd{d-l)"-'^ Igd 



n+l n 

A= Yl^/M'i2^' = yM^-2^"+i • JjyM^-2^^ 

i=l i=l 

<{d- l)"/2 . 2('^-l)" . rf"/2(d - l-jn{n-l)/2 . 2^^^-'^)' 

^ 2{n+l)d(d-l)"-'^+{n2+n)lgVd 
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Moreover, (#Q„+i) < 2d"+\ (#Q,) < 2(d - l)"+\ and so 

n+l n 

Q = n (#^^)''' = i*Qn+i?'-^' n (#^*)*'' 

n 

We apply (fTO]l using the previous inequalities, and get 

m ~ 
To assure that the minimum is attained inside the simplex, we apply a transformation that 
preserves the degree, but the bitsize of the polynomial is now bounded by r + 1 + dlg(n). 
Replacing this in the previous inequality, we get ^ < ■;^ — , where 

1 _ 2(n^+n)lgV^+(2+3n+d+(n2+3n+l)lgd _ 2{"+l)'^lg")'^('^-l)""^ . 2("+l)^'^('^-^)""^ (26) 

If we know that the system is zero dimensional then we could use Th. [3l Of course this is 
not always the case, hence we state the following bound, using ([5]), just as a reference. 

J_ < ^ ^ 2(("+l)^+"+^+("^+3"+l)lSf^)a!(c'-l)"~^ C27) 

m ~ rriDMM 

Let us compare the m-DMMp with other bounds that appear in the bibliography. In [21 Sec. 2, 
Rem. 2.17], the following estimation was computed, 

1 _ 22"+^n(i"+H-r+8nrf)j^2"+5d"+2n i2"+5d"+in2 

"iBLR (28) 

_ 22"+^nTd"+i+2"+5nd"+i(2n(i+dlgn+nlgd) 

which also holds with no assumption, but it is looser than r7iDMM„- 



p 



In [5] the authors derive a bound for the minimum of the absolute value of a polynomial, 

m — msY ' 



^ < ^, i.e. 



— = ((n + 2)2e"+3)(«+i)("+2)d"+i (n"(7i+ 1)^20^"+')''". (29) 

TTlBY 

The authors use the terminology evaluation hound for their bound. It holds when there is a 
0-dimensional projection; they prove that this is always the case for 
In |16j the following bound was computed: 

1 1 



— < 

m mjp 



2(r+l)d"+^^(n+l)d"+i^ (30) 



which has no restriction on the corresponding polynomial system. It is comparable to mDHH^j in 
general, but strictly looser when d > n. 

Example 12. Let us compute a lower bound on the value of f = {x -\-2y — 3) + [x + 2y — 4) , 
where d G {2,8,32}. Tiie polynomial is positive as it is a sum of squares. Consider the ideal 
I = if -z,fx,fy) C Z[x,y,z]. If {(1,(2,(3) belongs to the zero set of If, then iCsl > 2"'', b > 0. 
In Tab. [I] we present the estimations oflgb by the previous bounds. The true value is 6 = 0. 
When the degree is comparable to the number of variables (d = 2), then our bound and rrijp 
are comparable. When d > n, e.g. d = A and d = 32, then m^MMp is better than mjp by an order 
of magnitude. 

12 



5 Subdivision algorithms 

We use our results to bound the number of steps that any subdivision algorithm performs to 
isolate the real roots of a well-defined polynomial system. Then, we bound the complexity of 
Milne's algorithm in 2d. Our analysis can easily be extended to IR", however it is not clear 
what is the exact bit complexity of the elimination steps needed. 

We use DMM„, Th. [3l and Eq. ([1|) & (l9|), to bound the number of steps of a subdivision 
algorithm to isolate the real roots of a well-defined polynomial system as in ([T|). We assume 
the existence of an oracle that counts the number of real roots of the system inside a box in 
Q". Our aim is to compute the number of calls to the oracle in order to compute isolating 
(hyper-)boxes for all real roots. Realizations of such oracles are in [221 \2E[ 124] . see also jj3j. 

Suppose all roots of the system lie in a hypercube of side C, see Th. [3l At step h of the 
algorithm, the oracle counts the number of roots in hypercubes of side C/2^. We consider the 
whole subdivision algorithm as a 2"— ary tree T, where at each node we associate a hypercube, 
and to the root of the tree we associate the initial hypercube. Let #(T) denote the number of 
nodes. We will prune some leaves of T to obtain tree T' where it is easier to count its nodes. 

We proceed as follows. If t; is a leaf and has a sibling that it is not a leaf, then we prune 
V. If ui, . . . ,Uk, for some positive integer k, are leaves and siblings, such that they have no 
sibling that is not a leave, then we prune all of them except one that possess a hypercube that 
contains a real root. Notice that there is always at least one such node in ui, . . . ,Uk, because 
otherwise, the subdivision process in this path would have stopped one level before. If there 
exists more than one such node in ui, . . . ,Uk, then we keep arbitrarily one of them. It holds 
that #(r) < 2"#(r'), and we wiU count the nodes in T' . 

Each leaf of the tree contains contains a hypercube that isolates a real root of the system, and 
if there are at most R real roots, this also bounds the number of the leaves of T' . The hypercubes 
that correspond to the leaves of the tree have diagonals that are at least Aj = \^j — 7c |, where 
7c. is the root closest to jj. The length of their edges is at least \^j^i — 7c ,«|, where 1 < i < n. 
It holds that Aj = |7j — 7c.| > \jj^i — 7c ,i|, for any index i. The number of nodes from a leaf 



to the root of the tree is 



log^ 



Hence the number of nodes in #{T') is 



-ft r ^ -] -ft 



J = l 



#{T') = Y, log^ <R + RlgC-lgllAj. (31) 



j=i 



To bound the various quantities that appear, we will rely on Eq. Q and Th. [3l If the total 
degree of the polynomials is bounded by d, and ||/i||oo < 2"^, then IgC < nTd'^~^. To bound 
f3 • j^ Aj we use Eq. (j3]) of Th. [3] with £ = R. The hypotheses of the theorem, concerning the 
indices of the roots, are not fulfilled when symmetric products occur. In this case, we factorize 
quantity as ni=i ^i — Y[i=i '^i Y[i=i ^«) where Ri + R2 = R and the factors are such that no 
symmetric products occur. Then 

nf=iA. = nf=\A.nEiA.> 

2-R-{D-l){D+2) tf^ (j\2-2D-R j^-{n-l){2D'^+D{R+2)+R) _ 

If we take into account that R < D < d^, then 

-lognj=i Ai < 20"^ + 3D\gC + 3Dlgh + 5nD^lgB 

< 8(lg n + nlg d)d^'^ + 3n(n \gd + r)(i2"-i , 
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and for the total number of nodes of T' we have 

(#T0< i?+i?igc-ignf=iA, <z?+I)lgc-lg^f=lA, 
< 2d"(nTd'^-i) + S{\gn + n\gd)(f'' + 3n(nlgd + r)^^"-! 
= 0{n{n + d + T)(P''-^), 

and hence (#r) = 0(2" n(n + d + r) d^^^i). 

Theorem 13. Consider the polynomial system formed hy the polynomials in (T^. The number 
of steps that a subdivision algorithm performs in order to compute isolating boxes for all the 
real roots of the system is 0(2" D{D + lg C)) or 0(2" {d + r) d^""!). 

Remark 14. If we specialize n = 1 in the previous theorem, then we deduce that the number 
of steps of subdivisions algorithms for real root isolation of univariate integer, not necessarily 
square-free, polynomials is 0{d'^ Igd + dr). The optimal bound is 0{d'^ + dr) jTTJ. 

We now bound the complexity of Milne's algorithm [22] for isolating all real roots of a 
bivariate polynomial system. Milne's, so-called, volume function realizes the required oracle, 
see [15t I32j for experimental results. By SR(/, 5) we denote the signed polynomial remainder 
sequence of /, g. 

Proposition 15. jlBl [I^ We compute SQ{f,g), any polynomial in SIl{f,g), and Res{f,g) 
wrt X in Osiqip + q)^^^dT). The degree ofSIl{f, g) in yi, . . . ,yk is 0{d{p + q)) and the bitsize 
is 0{{p + (?)t). We can evaluate SR(/, 5) at x = a, where a G Q U {00} and £ (a) = a, in 
deiqip + q)^^^dmaii{T,a}). 

Let f,g € l\x,y\ with total degrees bounded by d and bitsize bounded by r. We are 
interested in isolating the real roots of the polynomial system f{x,y) = g{x,y) = 0, which we 
assume to be 0-dimensional. We introduce new parameters u, a, b and we eliminate a, b from 
the polynomials {/(a, b),g(a, b),V = u + [x — a){y — 6)}, where V is the volume function. After 
elimination, we obtain a polynomial h G (Z[a;, y])[ii]. We compute the Sturm sequence of h and 
its derivative w.r.t. u, hu, and we evaluate the sequence over u = 0. We obtain a sequence of 
bivariate polynomials in x, y. Now consider a box in the plane. We evaluate the sequence on 
each vertex of the box, and we count the number of sign variations. The number of real roots 
inside the box is j the sum of the sign variations |22j . 

We perform the elimination using iterated resultants. Using Prop. [15] we compute hi = 
ReSa{f{a,b),V{u,x,y,a,b)) G X[u,x,y,a,b] in Osid'^T). The total degree of hi is 0{d'^) 
and C{hi) = Oi^dr). Similarly, we obtain polynomial /i2 = ReSa{g{a,b),V{u,x,y,a,b)) G 
Z[u,x,y,a,b]. Finally, h = Resfe(/ii, /12) G Z[x,y,u] is computed in Osid^'^T). The degree of h 
in u is 0{d^) since the resultant of hi, /i2 has the factor n'^g(/(^'0))^g(9(^'0)) = u'^ . The degree of 
hmx,y is dsid^) and £ (/i) = 6{d^T). 

We compute the signed polynomial remainder sequence of h, hu and evaluate it at 0. This 
costs OB{d^^T). The evaluated sequence contains 0{d^) polynomials in Z[x, y] of degrees Osid^) 
and bitsize ©^(d^r). Each polynomial in the sequence is evaluated over a rational number of 
bitsize a in ©^(^^^(t + da)), and thus all of them in ©^(^^^(t + da)). 

In the worst case, a equals the bitsize of the separation bound, i.e. 0{d^T). Hence, the 
evaluation of the sequence costs Osid'^^T). Th. [13] indicates that we need to perform this 
evaluation 0{d'^ Ig d + d'^r) times. 

Theorem 16. Let f,g£ 'L[x,y] with total degrees bounded by d and bitsize bounded by r. 
Using the algorithm of Milne \22] . we can isolate the real roots of the system f = g = in 

OBid^'^T + d^^T^). 
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bound 


(d,r) = (2,5) 


(8,20) 


(32,85) 


12], Eq. dlHD 


1 lg("^BLR) 


27136 


6 684672 


1604321280 


[5], Eq. (I29D 


|lg("JBY)| 


1192 


74 000 


4696811 


[IS],Eq. (30) 


|lg("Jjp)l 


72 


15 360 


3 309 568 


Eq.dMD 


|lg(?TlDMMp)| 


87 


7457 


442 447 


Eq.dSZD 


|lg(mDMM)| 


54 


5 201 


324 506 



Table 1. Comparison of (the bitsize of) various bounds on the minimum value of the polynomial f — (x + 2y — 
3)'' + {x + 2y - 4)'', for d G {2, 8, 32} and r G {8, 20, 85}, resp. The bounds hold for all polynomials with same 
characteristics. 

Bounds on mutli-point evaluation of multivariate polynomials [23] could save at least two 
factors in the previous theorem. 
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